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Search for microwave emission from ultrahigh energy cosmic rays 
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We present a search for microwave emission from air showers induced by ultrahigh energy cosmic 
rays with the microwave detection of air showers experiment. No events were found, ruling out a 
wide range of power flux and coherence of the putative emission, including those suggested by recent 
laboratory measurements. 



I. INTRODUCTION 

Evidence for radio emission in the GHz regime from 
test-beam induced air showers by Gorham et al. [l| sug- 
gests a novel technique to detect extensive air showers 
(EAS) produced in the atmosphere by ultra-high energy 
cosmic rays (UHECRs). 

Microwave bremsstrahlung emission [l[ is expected to 
be isotropic and follow the EAS longitudinal develop- 
ment, providing an analog to the already successful flu- 
orescence detection technique @, Q , but with a tenfold 
increase in detector duty cycle and the elimination of at- 
mospheric systematics. In addition, thanks to the com- 
mercialization of the GHz frequency band, off-the-shelf 
components are readily available, allowing for inexpen- 
sive instrumentation of very large areas. A new detec- 
tor of this kind could allow current experiments like the 
Pierre Auger Observatory [|| and Telescope Array Q , as 
well as a future larger scale EAS observatory, to collect a 
statistically significant sample of UHECR measurements 
sensitive to composition at the highest energies which is 
crucial to unveil the origin of UHECRs @ . 

Several complementary approaches to microwave de- 
tection of EAS are currently being pursued, including 
the AMBER, EASIER and MIDAS detectors at the 
Pierre Auger Observatory, and the CROME experiment 
at KASCADE Also, new laboratory measurements 
are being performed || [Io| • Prior to its installation at 
the Pierre Auger Observatory, we have built and oper- 
ated at the University of Chicago campus a prototype of 
the Microwave Detection of Air Showers (MIDAS) de- 
tector |ll|, which is sensitive enough to detect EAS mi- 



crowave signals as estimated from measurements in [l[ 
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In this Letter, we report results of this first period of 
data taking and the corresponding limits on the power 
flux and coherence of microwave emission from EAS. 



II. MIDAS DETECTOR 



The MIDAS detector consists of a 4.5 m diameter mo- 
torized parabolic reflector deployed at the University of 
Chicago, instrumented with a 53-pixel prime-focus re- 
ceiver array built from commercially sourced low noise 
block feeds (LNBFs, C-Band 3.4-4.2 GHz). The corre- 
sponding field of view is approximately 20° x 10°, with an 
average system temperature of 65 K. Each pixel's signal 
is passed through a band-pass filter which rejects most 
of the anthropogenic noise, fed into an analog power de- 
tector and then digitized by a 20 MHz analog-to-digital 
converter (ADC). Self-triggering of the MIDAS detector 
is implemented through field programmable gate arrays 
(FPGA) in digital electronics boards. For each pixel, 
a first level trigger (FLT) is issued when a 1 /is ADC 
running sum over the pixel's data time stream is over 
threshold. The pixels' thresholds are dynamically regu- 
lated to keep the FLT rate close to 100 Hz. The second 
level trigger (SLT) searches for 4-pixel patterns of FLTs 
with a track-like topology across the camera and overlap- 
ping in time within the 10 fxs FLT time gate, as expected 
of EAS induced signals. When a SLT is found, a 100 
/is ADC time trace is written to disk for each pixel of 
the detector. SLT rates during data taking varied be- 
tween 0.01 Hz and 2 kHz depending on noise conditions. 
A data taking run lasts six hours, and is automatically 
restarted. Remote operation and monitoring of the MI- 
DAS detector is implemented. Details can be found in 

nam. 
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TABLE I. Table of cuts used in search program and their effect on selected data sample 



Cut 


Events Remaining 


After Cut 


(1) 


Less than 3 FLT pixels outside the SLT time window 


625,012 


(2) 


All SLT patterns are time-ordered down-going 


4,112 


(3) 


SLT pattern crossing time greater than 400 ns 


1,432 


(4) 


Traces in triggered SLT patterns contain only 1 pulse > 5 a 


979 


(5) 


Pulses > 5 a have a shape consistent with power detector's time constant 


924 


(6) 


FLT pixels matching a 5-pixel pattern topology with down-going time order 


21 


(7) 


Visual inspection of candidate events 






III. MICROWAVE SIGNAL SIMULATION 

The sensitivity of the MIDAS detector to EAS has been 
estimated with an end-to-end Monte Carlo simulation 
which includes a full description of the antenna charac- 
teristics. Following we parameterize in the simulation 
the microwave flux of the EAS at the MIDAS detector, 
If, as: 

if = if.ref ■ (p/po) ■ (d/R) 2 ■ (N/N rci ) a (1) 

where If !rc f is the power flux at a distance d = 0.5 m 
from a reference shower of E re j = 3.36 x 10 17 eV as 
given by R is the distance between the detector and 
the EAS segment, p (po) is the atmospheric density at the 
altitude of the EAS segment (at sea level), and N is the 
number of shower particles in the EAS segment. iV re f is 
the average number of shower particles at the maximum 
of the EAS development for a proton primary of energy 
E re f. N and N le { are given by a Gaisser-Hillas [lU pa- 
rameterization of their respective EAS. The parameter a 
characterizes phenomenologically the coherence scaling 
relationship for the EAS microwave emission. 
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FIG. 1. Simulated signal from an UHECR shower of energy 
3 x 10 18 eV at 5.5 km distance from the MIDAS detector with 
/f,rcf = 1-85 x 10~ 15 W/m 2 /Hz and quadratic scaling, a = 2. 
The histogram shows the FLT running sum calculated over 
the time trace for the triggered pixels of the camera (shown 
in the insert), with the horizontal line indicating the FLT 
threshold level. 

The microwave flux received at the MIDAS detector 
is propagated through the geometry of the reflector and 



LNBFs according to a microwave optics simulation, and 
converted to an ADC value. The gain pattern of the 
antenna used in the simulation was cross-checked with 
measurements of astrophysical GHz sources passing in 
the field of view. Absolute calibration of the MIDAS 
sensitivity with solar flux measurements was used for the 
ADC conversion. The simulated ADC signal is then pro- 
cessed using software triggering with the same logic im- 
plemented in the MIDAS FPGA [ll|. 

Coherent plasma emission can occur by the superposi- 
tion of emission wavefronts caused by the spatially coher- 
ent distribution of individual emitters. In this work, we 
used a phenomenological parameterization of the plasma 
coherence which may account for the different conditions 
present in an EAS and in the laboratory measurements. 
We assume a plasma of electrons to be composed of 
M subgroups of /i e coherent emitters, N c = M/i e - The 
total emitted power is then given by the incoherent sum 
of the emission from M subgroups: 

/f,tot = M/4r f ,i, (2) 

where if i is the power flux from a single electron in the 
plasma. Equation [5] can be written as 

if, tot = N°I ttl , (3) 

with the scaling parameter a = 1 + 1°^^ ■ The number 
of coherent emitters, p , and thus the value of a, will 
depend on the coherent spatial properties of the plasma, 
including its temperature, density, and electron energy 
distribution, which may significantly differ between an 
EAS induced plasma and a plasma generated with test- 
beams in laboratory measurements. In any case, a power 
law dependence as in Eq. [TJ is expected, from Eq. [3] and 
the proportionality of the number of plasma electrons 
to the number of shower particles. Since there is little 
understanding, both theoretically and experimentally, of 
the properties of the tenuous air plasma generated by 
EAS, we have no definite prediction for a, which may 
take any value from one (incoherent emission) to two 
(fully coherent emission). When comparing simulations 
to data, the degeneracy in Eq. [1] between If , rc f and the 
coherence scaling parameter a will be taken into account. 

Laboratory measurements [l[ suggest a reference power 
flux 7 f ° ref = 1.85 x 1CT 15 W/m 2 /Hz and full coherence at 
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shower maximum, a — 2. With these values of emission 
parameters, MIDAS should detect several tens of highly 
distinctive events per month. An example of a simulated 
event is shown in Fig. [TJ 



IV. DATA SAMPLE AND EVENT SELECTION 

The MIDAS detector took data during several months 
in 2011, under various conditions of anthropogenic noise 
and with several periods dedicated to calibration mea- 
surements. To ensure quiet and stable data taking con- 
ditions, we limited the search for EAS candidates to data 
runs with an average SLT rate less than 0.7 Hz. This 
requirement eliminates data taking periods with partic- 
ularly noisy conditions. The final data sample contains 
1.1 x 10 6 SLT events collected during 61 days of live time. 

We expect only sa 1600 SLT events from accidental 
coincidences of background fluctuations, indicating the 
data sample is likely dominated by anthropogenic noise. 
Simple event selection criteria, based on timing and trig- 
ger information expected from EAS geometries and labo- 
ratory measurements of detector electronics, were found 
to eliminate most of the noise events. Table U summarizes 
the selection criteria, and the corresponding number of 
selected events. 

With cut (1), we require that no more than two FLT 
pixels in the event not participating in any SLT pattern 
are triggered outside a time window defined by the small- 
est and the largest FLT time amongst the pixels forming 
SLTs. This cut rejects a class of background events char- 
acterized by a large number of camera pixels having a 
FLT, but distributed in time outside the SLT window 
and having topologies that do not form a SLT pattern. 

A large background rejection is obtained by asking that 
the time ordering of the FLT pixels is consistent with 
a down-going track topology for triggered SLTs, as ex- 
pected for EAS (cut (2)). Also, we require the FLT 
time difference between the latest and the earliest pixel 
of the SLT pattern to be larger than 400 ns (cut (3)), 
in order to reject a class of background events with co- 
incident pulses within the typical response time of the 
power detector. Cut (3) will also eliminate cosmic ray 
events which produce fast signals, like those with geome- 
try pointing towards the telescope, which would be in any 
case very difficult to distinguish from the overwhelming 
background of fast anthropogenic transients. In addition, 
each pixel participating in the SLT must have recorded 
only one pulse > 5 a (cut (4)), to reject events with mul- 
tiple pulses in the time trace, typical of noise induced by 
airplane altimeters and communications. With cut (5), 
we ensure that the detected pulse in each SLT pixel time 
trace is consistent with an RF pulse passing through the 
power detector - a pulse with a fall time shorter than the 
33 ns time constant of the power detector is likely to be 
a random fluctuation or anthropogenic noise. Finally, we 
search the events for a 5-pixel track topology with down- 
going time ordering (cut (6)). We expect 0.013 events 
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FIG. 2. One of the candidate 4-pixel events found in the 
data exhibiting topology which is consistent with expectations 
from an EAS. 

from accidental coincidences with such 5-pixcl topology 
in 61 days of live time, which would make any candidate 
event highly significant. Only 21 events survive this cut. 

A visual inspection of the selected events was then per- 
formed (cut (7)). All were identified as clear anthro- 
pogenic noise events, illuminating large portions of the 
camera with multiple coincident pulses that fall below the 
5 a threshold in cut (4). Background events with simi- 
lar topology and a stronger signals are seen frequently in 
the data and are indeed rejected by our selection criteria. 
The most likely explanation is that the residual back- 
ground events after cut (6) are produced by the same 
source of microwave radiation, but with a different ge- 
ometry relative to the detector. Notice that almost all 
of the 4-pixel events remaining after cut (5) are back- 
ground, as evident after a visual scan of the candidates. 
One example of a 4-pixel event compatible with an EAS 
signal is shown in Fig. [5] However, the characteristics 
of these events, which have signals close to the trigger 
threshold, are also compatible with a tail of the over- 
whelming background noise. No strong conclusion can 
be drawn on the origin of these kinds of events until a 
validation by a coincident detection at the Pierre Auger 
Observatory is performed. Thus, we use the stronger 5- 
pixel selection, which yielded a null result, to establish a 
limit on microwave emission from EAS. Notice that the 
large field of view of the MIDAS telescope and our search 
criteria, based only on topology and timing of the signals 
and with no requirement on the position of the maximum 
development of the shower, makes the result of our search 
insensitive to the composition of UHECRs. 



V. RESULTS 

In order to establish a limit on the microwave emis- 
sion from UHECR induced air showers we simulated a 
large sample of Monte Carlo events spanning a range of 
coherence scaling parameter a values between one and 
two, and of reference power flux Jf jre f values between 
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2.3 xlCT 16 W/m 2 /Hz to 4.6 x 1CT 15 W/m 2 /Hz (cf. Eq.Q} 
for the MIDAS configuration operating at the University 
of Chicago campus. For each pair of {a, if ilc f }, events 
are simulated in discrete bins of primary energies be- 
tween log 10 E = 17.65 and log 10 E = 20.05 in logarithmic 
steps of 0.1, sampled from an isotropic distribution. For 
each energy bin, the simulation is run until 4000 SLTs 
or 4 x 10 6 simulated EAS are reached, whichever occurs 
first. The result is then weighted with the UHECR en- 
ergy spectrum of [3] to produce a realistic energy spec- 
trum of simulated SLT events. 

The selection criteria (1) to (6) were applied to the 
simulated samples, and corresponding number of ex- 
pected events were derived. For example, for a fully 
coherent emission (a = 2) and a reference power flux 
-Tf.ref = If re f, we expect, in our 61 days of livetime, 
greater than 15 EAS candidate events above 2 x 10 18 eV 
passing all our search criteria. This is excluded at a 
greater than 5 a level by our null detection result. To 
cover the full space of {a,if ire f}, we used a linear sur- 
face interpolation of the grid of simulated samples. The 
shaded gray region in Fig.[3]is excluded with greater than 
95% confidence. We exclude the reference power flux 
i f ° ref as measured by [![, indicated by the solid horizon- 
tal line in Fig. [3l over a large range of partial coherence 
emission hypotheses. An incoherent emission for that 
reference power flux is not yet excluded by our measure- 
ment. 
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FIG. 3. Exclusion limits on the microwave emission from 
UHECRs, obtained with 61 days of live time measurements 
with the MIDAS detector. The power flux 7f irc f corresponds 
to a reference shower of 3.36 x 10 17 eV, and the a parame- 
ter characterizes the possible coherence of the emission. The 
shaded area is excluded with greater than 95% confidence. 
The horizontal line indicates the reference power flux sug- 
gested by laboratory measurements [l[. The projected 95% 
CL sensitivity after collection of one year of coincident oper- 
ation data of the MIDAS detector at the Pierre Auger Obser- 
vatory is represented by the dashed line. 



VI. CONCLUSIONS 

We reported the results of a search for microwave emis- 
sion from UHECRs with the MIDAS experiment. No 
candidate events were found in the 61 days of live time 
of the measurement, resulting in the strongest limits to 
date on the power flux and coherence of EAS microwave 
emission. In particular, we rule out fully quadratic scal- 
ing at the power flux value measured by Gorham et al. [l[ . 
This result has implications on the prospects for the mi- 
crowave detection of UHECRs, indicating a higher detec- 
tion energy threshold (or the need for a larger antenna) 
than the one so far considered. Notice that our limits 
apply to any microwave emission mechanism producing 
isotropic radiation that can be detected from the side of 
EAS at large distances, as in the fluorescence detection 
technique. The time response of the MIDAS electronics 
and the event search strategy is not optimized for detec- 
tion of very fast signals that may come from a beamed 
emission along the shower axis, as potentially observed 
in preliminary reports of EASIER Q and CROME 0. 

Measurements will continue, with the MIDAS detec- 
tor being deployed at the Pierre Auger Observatory in 
Malargiie, Argentina. Operating the MIDAS experiment 
in coincidence with the Auger Observatory will provide 
a much greater sensitivity to EAS events over the search 
reported here, eliminating the need for severe topolog- 
ical and timing cuts on candidate events. In addition, 



the rural environment of the site is significantly more ra- 
dio quiet than the urban environment of the University 
of Chicago campus, increasing overall detector livetime. 
The projected 95% CL sensitivity after one year of data 
is collected at the Pierre Auger Observatory is shown as 
a dashed line in Fig. |3l indicating that exclusion of inco- 
herent emission over a large range of reference fluxes will 
be achieved. 
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